Spin-isospin resonances: A self-consistent covariant description 
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For the first time a fully self-consistent charge-exchange relativistic RPA based on the relativistic 
Hartree-Fock (RHF) approach is established. The self-consistency is verified by the so-called isobaric 
analog state (IAS) check. The excitation properties and the non-energy weighted sum rules of two 
important charge-exchange excitation modes, the Gamow- Teller resonance (GTR) and the spin- 
dipole resonance (SDR), are well reproduced in the doubly magic nuclei '**Ca, ^"Zr and ^"^Pb 
without readjustment of the particle-hole residual interaction, 
exchange diagrams is demonstrated. 
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PACS numbers: 24.30.Cz, 21.60.Jz, 24.10.Jv, 25.40.Kv 



At present, spin-isospin resonances become one of the 
central topics in nuclear physics and astrophysics. Basi- 
cally, a systematic pattern of the energy and collectivity 
of these resonances could provide direct information on 
the spin and isospin properties of the in-mcdium nuclear 
interaction, and the equation of state of asymmetric nu- 
clear matter. Furthermore, a basic and critical quantity 
in nuclear structure, neutron skin thickness, can be de- 
termined indirectly by the sum rule of spin-dipole res- 
onances (SDR) [J Q or the excitation energy spacing 
between isobaric analog states (IAS) and Gamow- Teller 
resonances (GTR) Q. More generally, spin-isospin reso- 
nances allow one to attack other kinds of problems out- 
side the realm of nuclear structure, like the description 
of neutron star and supernova evolutions, the /?-decay of 
nuclei which lie on the r-process path of stellar nucleosyn- 
thesis [1, m , even the existence of exotic odd-odd nuclei 
and the efficiency of a solar neutrino detector 0]. 

It was realized long ago that the Random Phase Ap- 
proximation (RPA) is an appropriate microscopic ap- 
proach for charge-exchange giant resonances [1, H . The 
importance of full self-consistency was stressed [0] , and 
Skyrme-RPA calculations of charge- exchange modes ex- 
ist for about 30 years ■ Recently, a fully self-consistent 
charge-exchange Skyrmc-QRPA model has been devel- 
oped Self-consistency is an extremely important 
requirement for the analysis of long isotopic chains ex- 
tending towards the drip lines. On the relativistic side, 
so far the charge- exchange (Q)RPA model based on the 
relativistic mean field (RMF) theory has been developed 

However, the self-consistency of the RMF-f-RPA is not 
completely fulfilled for the following reasons. First, the 
isovector pion plays an important role in the relativis- 
tic description of spin-isospin resonances. Because of the 
parity conservation this degree of freedom is absent in 
the ground-state description under the Hartree approxi- 
mation. Therefore, the pion is out of control in this best- 
fitting effective field theory. Second, to cancel the contact 
interaction coming from the pseudovector pion-nuclcon 



coupling, a zero-range counter-term is needed with the 
strength g' ~ 1/3 exactly |15j . However, in order to re- 



produce the excitation energies of the GTR, g' must be 
treated as an adjustable parameter in RMF+RPA model 
with the value g' « 0.6 [l^.lTi|. 

In this Letter, for the first time a fully self- 
consistent charge-exchange relativistic RPA model is es- 
tablished, based on the relativistic Hartree-Fock (RHF) 
approach [l^, [l^. The two major advantages of this 
RHF+RPA approach are that the pion is included in 
both the ground-state description and the particle-hole 
(p-h) residual interaction, and that the zero-range pi- 
cnic counter-term with g' = 1/3 is maintained self- 
consistently. Without any adjusted p-h residual in- 
teraction or re- fitting process, we expect the present 
RHF+RPA approach to be reliable and to have predic- 
tive power. 

For a self-consistent calculation, the RPA p-h residual 
interaction must be derived from the same Lagrangian as 
ground-state [l^,[lB|- In the one- meson exchange picture, 
the interactions are generated by a-, uj-, p-, and tt- meson 
exchanges. The nucleon-nuclcon interactions read 

K(l,2) = -[g,7o]i[5^7o]2^a(l,2), (la) 
K.(l,2) = [5^7o71i[.9..7o7m]2^-(1,2), (lb) 
Vp{l,2) = [5p7o7^r]i-[5p7o7,,T]2Z?p(l,2), (Ic) 

K-(l,2) = -[— r7o757''i9fe]i • [— T7o757'5;]2 

xDAl,2), (Id) 

where 1)^(1,2) denotes the finite-range Yukawa type 
propagator 



A(l,2) 
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(2) 



As discussed before, a zero-range pionic counter-term 
with g' = 1/3 must be included, 

V;5(l,2) = .9'[— f7o757]i-[— ■?7o757]2'5(ri-r2). (3) 
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FIG. 1: (color online) Strength distribution of GTR in ^"^Pb 
calculated by RHF+RPA with PKOl (solid hue). The unper- 
turbed (HF) strength (dotted line), the calculation with only 
a + LO p-h residual interaction (dashed line) , and the calcula- 
tion excluding pion (/,r = 0) in the p-h residual interaction 
(dash-dotted line) are also shown. A Lorentzian smearing 
parameter F = 1 MeV is used. 



In the present RHF+RPA framework, both direct and ex- 
change terms must be taken into account and therefore 
the isoscalar mesons also play their role in spin-isospin 
resonances via the exchange terms. This is another dis- 
tinct difference from the RMF-I-RPA model. 

In order to verify the model self-consistency, we per- 
form a so-called IAS check according to the following 
property about IAS degeneracy: it is expected that 
the IAS would be degenerate with its isobaric multiplet 
partners if the nuclear Hamiltonian commutes with the 
isospin lowering T_ and raising T_|_ operators, which is 
true when the Coulomb field is switched off. While this 
degeneracy is broken by the mean field approximation, 
it can be restored by the self-consistent RPA calculation 
. As an example we calculate the IAS in ^"^Pb with the 
parametrization PKOl [l^ and we find the unperturbed 
excitations between —10.46 MeV and —8.96 MeV when 
the Coulomb interaction is put to zero, thus showing the 
isospin symmetry breaking. Then, the RHF+RPA cal- 
culation leads to i^iAS = 4 keV within the single-particle 
energy truncation [—A/, M+80 MeV] . This restoration of 
the IAS degeneracy indicates that the present approach 
is fully self-consistent. Furthermore, it should be empha- 
sized that the pion also plays its role in this restoration 
process. Therefore, the coefficient g' is not a free pa- 
rameter. Changing the value of g' would destroy the 
symmetry restoration process. For example, g' = leads 
to £;iAS = -801 keV. 

Taking the doubly magic nuclei ^^Ca, ^"Zr and ^"®Pb 
as test cases, and using the Gamow- Teller (GT) op- 
erator F^'^ = ^iC{i)T±{i), the excitation energies 
and strengths calculated with the fully self-consistent 
RHF+RPA approach using the parametrizations PKOl, 
PK02, PK03 [13 are summarized in Table [H These 



three parametrizations correspond to different sets of 
coupling strengths and meson masses in Eqs. ([l])-(l3]). In 
PK02 the pion is not included (i.e. f^^ = 0), whereas 
PKOl and PK03 have different constraints on the den- 
sity dependence oi Jtt- A good agreement with empirical 
energies is obtained without any re-adjusted parameter. 
All calculated strengths correspond to the main peak and 
contain 60-70% of the Ikeda sum rule. 

We can understand the different physical mechanisms 
between the present RHF+RPA and other RMF+RPA 
approaches by the following analysis. On the one hand, 
it has been shown that the ttNN interaction and its zero- 
range counter-term (g' w 0.6) are the dominant ingredi- 
ents in p-h residual interaction for the GT mode in the 
case of RMF+RPA [H, H!]. On the other hand, in the 
present RHF+RPA calculations, three parametrizations 
PKOl, PK02 and PK03 lead to similar resuhs for the 
GTR excitation energies. It should be emphasized that 
the pion is not included in PK02. This hints to the fact 
that the pion interaction is not the only dominant in- 
gredient for the GT excitations in this framework. The 
GT strength distribution in ^°*Pb with PKOl is shown 
in Fig. [TJ It is compared with the calculation in which 
the pion is excluded (/^r = 0) in the p-h channel, the cal- 
culation including only a + lu p-h residual interactions, 
and the unperturbed (HF) case. One can conclude that 
the isoscalar a- and w-mesons play an essential role via 
the exchange terms, whereas the pion just stands on a 
marginal position in determining the GTR strength dis- 
tribution. 

The relativistic RPA is equivalent to the time- 
dependent RMF in the small amplitude limit if the p-h 
configuration space includes not only the pairs formed 
from the occupied and unoccupied Fermi states but also 
the pairs formed from the Dirac states and occupied 
Fermi states [l^. Based on this idea, a relativistic re- 
duction mechanism of the Gamow- Teller strength due to 
the effects of the Dirac sea states was pointed out [l^ . 
This kind of reduction mechanism appears in both nu- 
clear matter ^ and finite nuclei [11, [13. In Table HH 
we show the contributions to the Ikeda sum rule values 
coming from the Fermi (5^) and Dirac (5^) sectors of 
both T± channels. It is explicitly shown that the Ikeda 
sum rule 



_ ^GT 



3(N- Z) 



(4) 



can be 100% exhausted only when the effects of the 
Dirac sea are included. The reduction factors, 1 — {S^ — 
Sl)/{S_ - S+), of ^^Ca, 9°Zr, ^ospb by the present self- 
consistent approach indicate to which extent the antinu- 
cleon degrees of freedom play a role. 

It has been proposed that the neutron skin thickness 
could be extracted via the spin-dipole sum rulefH, 
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with the SD 

(T{i)]j=0,1.2T±{i). 



operator 
While 



experimental results in 
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TABLE I: GTR excitation e nerg ies in MeV and strength in percentage of the 3(A'' — Z) sum rule within the RHF+RPA 
framework. Experimental [13, [ll [13, [13, [111 and RMF+RPA [ll results are given for comparison. 
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60.5 




PK03 


10.42 


70.7 


15.71 


68.9 
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67.7 



RMF±RPA DD-MEl 10.28 72.5 15.81 71.0 19.19 70.6 



TABLE II: Ikeda sum rule values from Fermi (S^) and Dirac (5°) sectors. S-, S+ are the sum rule values of the T- and T+ 
channels, respectively. The reduction factor, 1 — (Sl — 5'+)/(S_ — 5'+), is given in the last column. 
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FIG. 2: (color online) Strength distributions in both T± chan- 
nels of the SDR in ^"Zr calculated by RHF±RPA with PKOl. 
The dash-dotted, dotted, dashed lines show the 0-, 1-, 2- con- 
tributions respectively, while the solid line shows their sum. 
A Lorentzian smearing parameter F = 2 MeV is used. The 
experimental data are shown as filled symbols 0]. 

both ^°Zr(p, n) [l^l and ^°Zr(n,p) (2^ channels are now 
available, the SDR has become a hot topic. 

Both T_ and T+ channels of the SD strength distri- 
butions are shown in Fig. O Here, the horizontal axis 
is the excitation energy measured from the ground-state 
of the parent nucleus ^°Zr. The RHF+RPA calculations 
reproduce the strength distributions up to the giant res- 
onance region in both channels without any quenching 
factor, even though they show a more pronounced struc- 



ture than the experimental spectra. Especially the dom- 
inant resonance structure centered at k, 27 MeV in 
the T_ channel is well reproduced. It is also found that 
the three components of the SDR follow the same energy 
hierarchy E{2~) < E{1~) < E{0~) as that found in the 
recent self-consistent Skyrme-RPA calculations [1^ [l^] ■ 
Since the Ip-lh configuration space of the RPA calcula- 
tions is restricted, the discrepancy beyond the giant res- 
onance region can be understood. Furthermore, in the 
first part of Table Hill the neutron and proton rms radii, 
the neutron skin thickness and the SD sum rule value of 
^°Zr are listed. They all agree with the data within the 
experimental accuracy. 

For the spin-dipole sum rule values, it is worth to 
note that a reduction due to the effects of the antinu- 
cleon degree of freedom in analogy with the GTR case is 
found. The present calculations show that this reduction 
is around 6.5% in ^°Zr. 

Since the neutron skin thickness of ^'^®Pb is important 
in many aspects of nuclear physics and astrophysics, it is 
worthwhile to investigate the SDR properties of ^"'^Pb. 
From Fig. [3] and the second part of Table [1111 the present 
approach predicts a dominant resonance structure at E — 
20 ~ 30 MeV in the T_ channel and a small bump at 
i? = 5 ~ 10 MeV in the r+ channel. Furthermore, a 
reduction factor of about 5.5% is predicted by the three 
different parametrizations. 

In conclusion, for the first time, a fully self-consistent 
charge-exchange relativistic RPA model based on the 
RHF approach is established. The IAS degeneracy bro- 
ken by the RHF approximation can be accurately re- 
stored in the present self-consistent RPA calculations. 
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TABLE III: SD sum rule values and neutron skin thickness of Zr, Pb in RHF+RPA approach. Neutron and proton rms 
radii and corresponding data from SD experiment 0| are given for comparison. S{rn , rp) stands for the RHS of Eq. ([5]). 
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FIG. 3: (color online) Same as Fig. [J] but for the nucleus 
2'«Pb. 

Compared with RMF±RPA, the isoscalar mesons (cr, cu) 
are found to play an essential role in spin-isospin res- 



onances via the exchange terms. The GTR excitation 
energies and their strengths can be reproduced in the 
present self-consistent RPA calculation while maintain- 
ing g' = 1/3 in the contact counter-term. The SDR 
strength distributions up to the giant resonance region in 
both channels are well reproduced in ^"Zr. Furthermore, 
a spin-dipole sum rule reduction mechanism due to the 
effects of the Dirac sea is found. The SD reduction factor 
for ^"^Zr is around 6.5%. Finally, the dominant structures 
of the SD strength distribution in ^"^Pb are shown and 
a strength reduction of about 5.5% is obtained. 
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